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ABSTRACT 

Absorption lines from water, methane and carbon monoxide are detected in the atmosphere of 
exoplanet HR8799 b. A medium-resolution spectrum presented here shows well-resolved and easily 
identified spectral features from all three molecules across the K band. The majority of the lines are 
produced by CO and H 2 O, but several lines clearly belong to CH 4 . Comparisons between these data 
and atmosphere models covering a range of temperatures and gravities yield log mole fractions of H 2 O 
between -3.09 and -3.91, CO between -3.30 and -3.72 and CH 4 between -5.06 and -5.85. More precise 
mole fractions are obtained for each temperature and gravity studied. A reanalysis of H-band data, 
previously obtained at similar spectral resolution, results in a nearly identical water abundance as 
determined from the AT-band spectrum. The methane abundance is shown to be sensitive to vertical 
mixing and indicates an eddy diffusion coefficient in the range of 10 ® to 10 ® cm^ s“^, comparable to 
mixing in the deep troposphere of Jupiter. The model comparisons also indicate a C/0 between ^ 
0.58 and 0.7, encompassing previous estimates for a second planet in the same system, HR8799 c. 
Super-stellar C/0 could indicate planet formation by core-accretion, however, the range of possible 
C/0 for these planets (and the star) is currently too large to comment strongly on planet formation. 
More precise values of the bulk properties (e.g., effective temperature and surface gravity) are needed 
for improved abundance estimates. 

Subject headings: planetary systems - stars: atmospheres - stars: low-mass, brown dwarfs 


1. INTRODUCTION 

The HR8799 planetary system remains unique as the 
only system to have mu ltiple, directly imaged, planets 
(jMarois et al.ll200^l20ir)ll . Four planets orbiting IIR8799 
have been monitored regularly since their discovery, pro¬ 
viding the astrometric data needed to estimate their 
orbital properties and masses. The current astromet¬ 
ric d ata indicate that all four planets are less than 13 
MT,„.(lMarois et aLl2010l:lF abrvck v fc Murrav-Clavl201d : 
iCurrie et al.ll201l7 Puevo et alir2014ll . 

The planets are frequently observed for the purposes of 
characterizing their atmospheric properties and, in par¬ 
ticular, their chemical compositions. The two outer most 
planets, b and c, have the most comprehensive wave¬ 
length coverage including near-infr ared (near-IR) spec¬ 
troscopy at low, R 50 — 100 dBowler eUaL 1201(11 
iBarman et al.l l 2011 al FO openheimer et al. l l2013ri and 
medium, R ~ 4000 dKononackv et al.ll2nT^ resolutions. 
The Gemini Planet Imager (GPI) has measured low reso- 
lution spectra of planets c and d (iMacintosh et al.ll20l4 


lingraham et al.l[2014fl . Spectroscopy of multiple plan¬ 
ets orbiting HR8799 allows direct comparisons of atmo¬ 
spheric compositions for a coeval set of planets formed 
from the same protoplanetary disk. The only other plan¬ 
etary system for which such a comparison is currently 
feasible is our own Solar System. 

The giant planets in our Solar System were likely born 
with bulk compositions determined by their initial loca¬ 
tion in the planetary disk and the specific solid-to-gas 
(Mgoiid/Mgas) accretion history they experienced. As 
these initially hot planets cooled with time, their at¬ 
mospheres experienced various levels of vertical mixing 
and condensation. The atmospheric composition may 
have been further altered by continued accretion of solid 
bodies or mixing with a deep metal rich core. Con¬ 
sequently, the present-day mole fractions for important 
trace molecules (H 2 O, CH 4 , CO 2 , CO, NH 3 , N 2 ) of gi¬ 
ant planets in our Solar System are the result of numer¬ 
ous chemical and physical processes. Inferring elemental 
abundances from molecular abundances requires a clear 
understanding of the atmospheric chemical and dynam- 
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Figure 1. Top and bottom panels compare the medium-resolution near-infrared spectra of HR8799 b (red, this work) and HR8799 c 
(blue; Konopacky et al. 2013), both normalized to the same flux at 2.15 fim. The lower panel includes a portion of the methane absorption 
reduced to the same resolution and sampling as the observed data and offset for comparison. Regions of strong methane absorption are 
indicated with vertical dotted lines. 


ical history. Even for Jupiter and Saturn only upper 
limits on carbon-to-oxygen ratio (C/0) have been mea¬ 
sured, a poignant r eminder of how ch a llenging such mea¬ 
surem ents can be (|Wong et al.l 12004 iVisscher fc Fegle5^ 
I2005D . Despite this complex connection between present- 
day composition and the composition at birth, it re¬ 
mains plausible that the present-day molecular abun¬ 
dances hold clues pertaining to the formation history of 
our giant planets. By observing young giant planets, like 
those orbiting HR8799, billions of years of atmospheric 
evolution that blurs the connection between atmospheric 
properties and the formation process is avoided. 

Presented below are new observations of HR8799 b 
that provide comparable wavelength coverage, spec¬ 
tral resolution and signal-to-noise (SNR) as similar ob¬ 
servations of HR8799 c that revealed individual re- 
solved water and carbo n monoxide absorption features 
(jKonopackv et al.ll2013L hereafter K13). For HR8799 b, 
the mole fractions of these two molecules as well as 
methane are determined using well-resolved spectral fea¬ 
tures. From these mole fractions, the coefficient of eddy 


diffusion (RTzz) and C/0 are estimated. 

2. OBSERVATIONS AND SPECTRUM EXTRACTION 

HR8799 b was observed in 2013 o n July 25, 26 and 27 
(UT) with the OSIRIS instrument (jLarkin et al.l [2(10^ 
in the K band using an identical instr ument configu¬ 
ration as previous work by this group (iBarman et al.l 
l2011al hereafter Bll) resulting in 5.2 hrs of on-target 
integration time under good observing conditions. The 
data were calibrated and rectified to p roduce three- 
dime nsional basic calibrated data cubes (jKrabbe et al.l 
l2004l . The data were initially binned to R ^ 100 
and speckles suppressed following Bll. As in Bll, at 
low-i?, the spectrum SNR is limited by residual corre¬ 
lated speckle noise. Despite the increase in signal from 
combining the 2009 to 2013 observations, the updated 
low-i? spectrum has mostly unchanged uncertainties and 
matches the earlier Bll spectrum. The AT-band spec¬ 
trum was absolute flux calibrated using the magni¬ 
tude (Mk^ = 14.15 ± 0.1) reported in Bll. 

At the full resolution and sampling provided by 
OSIRIS, errors in flux are correlated across several wave- 
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Figure 2. Top: Continuum-subtracted spectrum for HR8799 b. 
Template spectra for pure CH 4 (orange), H 2 O (blue) and CO 
(green) are also plotted and offset by an arbitrary amount. Bot¬ 
tom: Cross-correlation functions for HR8799 b and the template 
spectra plotted in the top panel. 

length channels and may alter the relative depth of nar¬ 
row absorption features. Subtracting the continuum, 
modeled as a smoothed version of the observed spec¬ 
trum, removes many residual speckle artifacts that vary 
smoothly with wavelength (K13). Inevitably, however, 
speckle artifacts that depend more strongly on wave¬ 
length (e.g., produced by optics near the focal plane) 
likely remain even in a continuum-subtracted spectrum. 
Despite potential lingering artifacts, this strategy worked 
well in K13 and is adopted here using the ensemble of 
data from 2009 through 2013, median-combined into a 
single spectrum. Given the wide angular separation of 
HR8799 b, residual speckle artifacts should be less severe 
than they were for HR8799 c. The RMS uncertainties 
in the continuum-subtracted spectrum closely matches 
the photon-noise limit. The low-i? spectra from various 
observing runs show no significant differences across this 
four year time frame and intrinsic variability would likely 
occur over timescales of days rather than years and im¬ 
pact broader wavelength ranges than the spectral lines 
of interest. 
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Figure 3. Same as Figure[2l but for H band. Only very weak CO 
lines are present in the H band and, therefore, were not searched 
for. The observed continuum-subtracted spectrum shows no corre¬ 
lation with the CH 4 template. 

are important when studying atmospheric abundances. 

The opacity of solid and liquid particles suspended 
in the atmosphere (clouds) are included using the pa¬ 
rameterized intermediate cloud model described in Bll. 
Briefly, the lower boundary of the clouds are determined 
by chemical equilibrium while the upper boundary is de¬ 
scribed by an exponential decay that begins at a spec¬ 
ified pressure. This outer pressure is a free parameter 
that establishes the cloud thickness, allowing for a range 
of models with cloud opacity between the high (DUSTY) 
and cloud-free (CON D) cases often used to bracket the 
importance of clouds (j Allard et al.ll200ill . The particle 
size distribution is a log-normal with mean size set to 5 
/im. Cloud thickness plays an important role in deter¬ 
mining the overall spectral shape and is included as a 
free parameter in the model fits discussed below (follow¬ 
ing Bll). 

Synthetic spectra were calculated with a wavelength 
sampling of O.OSA from 1.4 to 2.5 ^m. Each of these 
synthetic spectra was convolved with a Gaussian ker¬ 
nel with FWHM matching the OSIRIS spectral resolu¬ 
tion before interpolating onto the observed wavelength 
grid. These medium-resolution spectra were continuum- 
subtracted following the steps in K13. 


3. MODEL SPECTRA 

For this study, the grid of exoplanet atmo¬ 
spheres models described in B ll was updated to in¬ 
clude the methane linelist from I Yurchenko fc Tennyson! 
(I2014D supplemented with optical opacities from 
IKarkoschka fc Tomaskol (120101 1. The former update 
added approximately 10 billion transitions to the over¬ 
all molecular opacities and improves the accuracy of 
line strengths at high temperatures - across the K- 
band alone there are roughly a billion transitions. The 
methane abundance is low in non-equilibrium chemistry 
models appropriate for HR8799 b (Bll) and, conse¬ 
quently, there were only minor changes to the model 
atmosphere structures after this opacity update. Nev¬ 
ertheless, the ExoMol methane list is currently the most 
accurate and complete list available and the observations 
analyzed here have a resolution where accurate line data 


4. RESULTS 
4.1. HR8799 b versus c 

Figure[T]compares the R ^ 4000 spectrum of HR8799 b 
to that of HR8799 c. Both exoplanet spectra contain 
a similar set of absorption features, with most of the 
similarities at wavelengths of prominent H 2 O absorption; 
however, most lines are deeper in the spectrum of b. The 
CO (2,0) band-head at 2.3 ^m is detected now with 
medium spectral resolution but only marginally present 
at low-i? (Bll). This band-head appears slightly deeper 
and has a shallower slope in the spectrum of b than c, the 
result of stronger water absorption and additional CH 4 
lines for HR8799 b in this wavelength range. Three re¬ 
gions of CII 4 absorption are visually identifiable between 
2.15 and 2.4 /im (lower panel of Fig. [T]). The strongest 
methane lines are seen around 2.32 and 2.37 /rm. None 
of these features are present in the spectrum of c, consis- 
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Figure 4. Example synthetic spectra (red, blue, and green) for 
three solar-abundance models. The range of model-observation 
agreement shown here is representative of most models explored 
here. The largest discrepancies are between 3.5 and 4.5/rm. See 
text for details of the model and references for the observations. 

tent with the methane non-detection reported by K13. A 
few additional CH 4 lines are marginally identified around 
2.25 /tm, however, many of the weaker CH 4 lines overlap 
with those of H 2 O and CO, making visual identification 
difficult. 

4.2. Cross-Correlation 

As an additional confirmation that spectral features 
from all three molecules are present in the OSIRIS spec¬ 
trum of HR8799 b, a cross-correlation analysis was used 
following K13. In this case, a new CH 4 template was 
calculated using roughly a billion of the strongest trans i- 
tions across the K band (jYurchenko fc Tennvsonll2014D . 
Wate r and CO template s we re calculated usin g lists 
from iBarber et al.l (j200(ifl and iGoorvitcI] (jl994[ ). Fig¬ 
ure [2] compares the cross-correlation functions (CCFs) 
for each template. Peaks in the CCFs are found for 
all three molecules at identical velocities, centered on 
zero km/s. As was the case for HR8799 c, a peak with 
maximum near I is found for H 2 O, indicating that this 
molecule contributes most of the spectral lines. CO also 
shows the characteristic ringing pattern, produced by 
near-repeating patterns of CO lines, many separated by 
roughly 200 km/s in velocity. This same CCF pattern 
was found for CO in HR8799 c by KI3, with peak of 0.6 
compared to 0.4 here. A peak in the CCF for CH 4 is 
present, but not as prominent as found for H 2 O, mainly 
because much of the weak-line information is blended 
with other stronger lines at this resolution. 

iJ-band observations were previously obtained using 
the same telescope and instrument combination as used 
for A'-band and were originally analyzed at low resolution 
in Bll. These data have been reanalyzed and continuum- 
filtered at full resolution (also R 4000) following the 
same steps as K band. The average SNR of these data 
is about 3 times lower than the K band data and show a 
pattern of lines that are visually difficult to identify. Bll 
concluded that water is the dominant molecular opacity 
source in this wavelength range and this is confirmed here 


by a cross-correlation analysis. The CCFs for H 2 O and 
CII 4 templates across the H band are plotted in Figure 
[3l Correlating the data with an i7-band H 2 O template 
yields a strong peak in the CCF (Fig. |3|). No peak was 
detected for an i7-band CH 4 template (using the ExoMol 
list). 

The HR8799 b AT-band spectrum also correlates well 
with the HR8799 c spectrum, another indication that 
both planet spectra share similar spectral patterns of 
H 2 O and CO lines. The CCF was also recalculated for 
HR8799 c using the new CII 4 template and no peak was 
found, confirming the non-detection of CH 4 reported by 
K13. 

4.3. Model Comparisons 

For the atmospheric abundance study presented be¬ 
low, it is important that the underlying thermal struc¬ 
ture of the models, especially across most of the pho¬ 
tosphere, is reasonably correct. Atmospheres of giant 
planets are generally close to local thermodynamic equi¬ 
librium and, thus, a model matching the planet’s spectral 
energy distribution over a wide wavelength range should 
have a thermal profile that reasonably approximates the 
planet’s average thermal structure. A similar model com¬ 
parison as performed by Bll was repeated, but compar¬ 
ing updated models with the new CII 4 list to the new 
medium-resolution unfiltered spectrum discussed above 
and additional flux calibrated data spanning 1 to 5 ^m. 
Fo r wavelength s less than 1.8 /rm the z/Y-band flux 
of iCiirrie et al.l ( 120111). the low -A P1640 J-band spec- 
trum (lOopenheimer et al.l 120131 ) and the i7-band spec¬ 
trum from Bll were used. Photo metric data was used 
at wavelengths loriger than 3 Mm (iGalicher et al.l [20111 : 

iSkemer et al.l l20ia iCurrie et al.l 120141 ). The data used 

for the model fits are plotted in Figure IH 

The relative flux calibration of the near-IR spectra 
from OSIRIS and P1640 could impact the model com¬ 
parisons, as these spectra are calibrated using H and ATg 
photometry t hat have a range o f reported values and 
uncertainties (iMarqis et all 120081 iMetchev et al.l 120091 
lEsposito et al.ll2013f) . It is also possible that HR8799 b 
is variable as indica ted by studies of brown dwarfs 
(|Metchev et al.l 1201511 . Such potential issues are not 
accounted for here and the H and K flux calibration 
described in Bll is used. A small scaling was ap¬ 
plied so that the i7-band portion of the P1640 spectrum 
matches the i7-band OSIRIS spectrum, resulting in a 
slight change in the absolute fluxes of the P1640 J-band 
spectrum plotted in Figure 21 

Example spectra of the best matching solar abun¬ 
dance models are plotted in Figure IH for a range of 
Teff between 900 and HOOK and log( 5 ) between 3.0 and 
4.0. The best matching model has the same gravity as 
found previously by Bll (log((?) = 3.5) but cooler (Tes 
= 1000K). Across the near-IR, the model spectra do a 
reasonable job reproducing the observations. When fit¬ 
ting the observations (photometric and spectroscopic), 
all data are weighted equally. Between 3 and 5^m, how¬ 
ever, the model comparisons to available photometric 
data (from a variety of telescopes, epochs and image 
processing methods) show disagreements at L' and ~ 4 
/rm. Similar levels of disagreement at this wavelength 
are seen in other rnodel comparisons (iMarlev et al1l2012t 

ISkemer et al.ll20l3 : iCurrie et al.ll2014l) . These nhotomet- 
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ric observations probe similar atmospheric depths and, 
therefore, pull model fits in opposite directions in both 
Teff and log( 5 ). Further observations across this spec¬ 
tral region is probably warranted. Previous studies of 
HR8799 b result in a range of Tea between 750 and 1200K 
and log{g) between 3 and 4.5 (see Marley et al. 2012, 
their Table 1), encompassing the range explored here. 

Model atmosphere parameters are sometimes inten¬ 
tionally biased to match predictions of cooling tracks, 
providing guaranteed consistency between Tea, log(( 7 ), 
(and Lboi) and mass/radius expectations. Such forced 
agreement can be useful, but is not done here. More im¬ 
portant is finding a range of atmospheric temperature- 
pressure profiles that yield pseudo-continua that are as 
consistent as possible with the observed spectral energy 
distribution. Good agreement between bulk parameters 
found by atmosphere-only fitting and those inferred from 
cooling tracks, age and luminosity can be challenging 
for low-temperature planet-mass objects, likely a result 
of complexities of their atmospheres not reproducible 
by time-independent one-dimensional models with sim¬ 
plistic cloud prescriptions. Consequently, a conservative 
ranges of log(g) (± 0.5) and Teff (± lOOK) around the 
best matching values are used for the abundance analy¬ 
sis discussed below, despite potential inconsistencies with 
cooling track predictions at the high and low ends of 
these ranges. 

4.4. Mole Fractions of CO, CH^ and H 2 O 

At low spectral resolution, the slopes on either side of 
the AT-band peak are very similar to those seen in sub- 
stellar atmospheres where water is the dominant opacity 
source. Bll showed (their Fig. 14) that CO and CH 4 
likely contribute to the AT-band opacity but, at low res¬ 
olution, the CO band head was not confidently detected 
nor any significant hints of CH 4 , making it di ff icult to as¬ 
sess their relative importance. iBowler et al.l (|2010ll also 
found no evidence of methane absorption in their narrow- 
band spectrum covering ~ 2.1 to 2.2 /rm. With resolved 
absorption features from all three molecules identified, 
their individual abundances can now be determined. 

Evidence for quenching of CO and CH 4 in the at¬ 
mospheres of t he HR8799 planets has been estab¬ 
lished elsewhere (iHiiiz et al.ll201 0HBarnian et al.ll2011al : 
iSkemer et al.l 1201 2t iMarlev et al.l 1201^ An impor¬ 
tant consequence of this mixing-induced non-equilibrium 
chemistry is that the mole fractions should be nearly in¬ 
dependent of height across the photosphere. The absence 
of height-dependence greatly simplifies the parameteriza¬ 
tion of the mole fractions, allowing straightforward model 
fitting using standard minimization. 

The mole fractions are determined sequentially by 
computing grids of synthetic spectra (continuum fil¬ 
tered) with scaled molecular abundances for 9 different 
temperature-pressure (T-P) profiles corresponding to Teg 
= 900, 1000, and HOOK and log(g) = 3.0, 3.5, and 4.0. 
These temperatures and gravities bracket those found 
when comparing models to the broad wavelength data 
shown in Figure H) A spectrum is calculated for each 
mole fraction value and model temperature-pressure pro¬ 
file. The mole fractions of CO, CH 4 and H 2 O are each 
scaled from 0 to 1000 relative to their non-equilibrium 
solar abundance values using a uniform logarithmic sam- 



log,o Mole Fractions 


Figure 5. distributions for H 2 O, CO and CH 4 mole fractions 
when fitting model spectra with refj= lOOOK and gravity equal 
10® cm sec“^ to the K (solid) and H (dashed) bands. Vertical 
dotted lines indicate the best-fit mole fractions and vertical dashed 
lines indicate solar (non-equilibrium) mole fractions. CH 4 was not 
detected in the ff-band data, with a 3—cr upper limit of 10~® ®. 

pling (resulting in 100 synthetic spectra per T-P profile). 

The mole fraction of H 2 O was determined first, keeping 
CO and CH 4 at their solar values. To avoid biasing in 
the H 2 O mole fractions caused by overlapping CO or 
CH 4 lines, the fit was restricted to wavelengths less than 
the CO band head, while masking the strongest CH 4 
line region (near 2.2 p,m). Following this step, a second 
grid of synthetic spectra was calculated with the H 2 O 
mole fraction equal to the best-fit value from the previous 
step, CH 4 at its solar value, and now mole fractions of 
CO scaled as previously done for H 2 O. In this case, only 
wavelengths greater than (and including) the CO band 
head were fit, while masking the strong CH 4 line near 
2.32 and 2.37p,m. In the final step, using the best-fitting 
H 2 O and CO mole fractions, a third grid of synthetic 
spectra was made with scaled CII 4 mole fractions. Only 
the strongest CH 4 lines from 2.2 to 2.37 /im were included 
in this fit. 

The model with the overall lowest from this three- 
step fitting process has a slightly lower gravity (log(g)= 
3.0) than the best matching solar abundance model com¬ 
pared to the full SED. The x^ distributions for the mole 
fractions in this model are plotted in Figure [SJ with log 
mole fractions of H 2 O, CO and CH 4 are —3.85 ± 0.01, 
—3.67 ± 0.02, and —5.85 ± 0.04, respectively, with the 
model spectrum compared to the AT-band observations 
in Figure [51 Both H 2 O and CO are about half the 
non-equilibrium solar value while CH 4 is about 15 times 
lower than the non-equilibrium solar value (see vertical 
dashed lines in Fig. [S]). The best matching mole frac¬ 
tions for each T^g and log(g) are listed in Table 1, with 
formal 1 -cr errors determined from the x^ distributions. 
These values can be compared to the range of solar (non¬ 
equilibrium) mole fractions by using Figure [12] (see ver¬ 
tical dashed line), discussed below. 

The filtered A'-band spectrum is used to determine 
the molecular abundances primarily because it is a single 
data set from a single instrument with uniformly char- 
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Figure 6. Top: Continuum-subtracted spectrum for HR8799 b (black) compared to refj= lOOOK and log{g)= 3.0 model with abundances 
equal to those found in Figure |5] Bottom: Spectrum of HR8799 b (black) compared to same model plotted above but with continuum 
intact. 


acterized uncertainties. It is important, however, to ver¬ 
ify that the resulting mole fractions still yield spectra 
that are in reasonable agreement with the full SED. Fig¬ 
ure [7] compares model spectra to the full SED for three 
different and log(g) values and corresponding mole 
fractions listed in Table 1, and the comparisons remain 
good across the near-IR. The thermal IR includes ab¬ 
sorption bands from CH 4 , CO and CO 2 , providing a po¬ 
tential secondary test of the abundances. In all cases, the 
model flux at 3.3/rm, which probes a methane fundamen¬ 
tal absorption band, is close to the observed value. The 
fluxes across the CO and CO 2 absorption bands between 
4 and 5/im changes considerably for the range of values 
found here, however, the observed ground-based photom¬ 
etry, with measurements on either side of the absorption 
bands, do not probe CO or CO 2 absorption well. 

As discussed above, H 2 O is detected in the med- 
resolution if-band spectrum by cross-correlation (Fig. 
[SI). Mole fractions for H 2 O were obtain by fitting the 


iJ-band data in a similar manner as for K band, using 
synthetic spectra with scaled H 2 O and CH 4 mole frac¬ 
tions. A clear x^-minimum was found for water for the 
same mole fraction as found by fitting the AT-band spec¬ 
trum, but with much larger uncertainties. CH 4 was not 
detected (consistent with the cross-correlation test) and 
only an upper limit was found (see dashed lines in Fig. 

El). 

The mole fractions determined above are sensitive to 
the thermal profile, surface gravity, individual line broad¬ 
ening parameters, and potential residual artifacts in the 
data. Uncertainties associated with these properties are 
difficult to quantify without a more complex statistical 
analysis of the model comparisons and is beyond the in¬ 
tended scope of this paper. For example, the mole frac¬ 
tions depend on the model atmosphere surface gravity, 
but it is unlikely the planet’s gravity falls outside the 
range explored here, given the limits placed on grav¬ 
ity by other data. The details of the cloud properties 
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Figure 7. Same observations as plotted in Figure[4]but compared 
to model spectra with non-solar H 2 O, CO and CH 4 mole fractions 
(see Table 0. 


Table 1 

log Mole Fractions 


molecule 

Tea 

3.0 

log(ff) 

3.5 

4.0 


900 

-3.82±0.01 

-3.39±0.01 

-3.09±0.01 

H 20 

1000 

-3.85±0.01 

-3.55±0.01 

-3.12±0.01 


1100 

-3.91±0.01 

-3.58±0.01 

-3.24±0.01 


900 

-3.61±0.02 

-3.48±0.02 

-3.30±0.02 

CO 

1000 

-3.67±0.02 

-3.61±0.02 

-3.33±0.02 


1100 

-3.72±0.02 

-3.61±0.02 

-3.30±0.02 


900 

-5.85±0.04 

-5.55±0.04 

-5.18±0.04 

CH 4 

1000 

-5.85±0.04 

-5.48±0.04 

-5.12±0.04 


1100 

-5.79±0.04 

-5.42±0.04 

-5.06±0.04 


(coverage and thickness) are partially mitigated by the 
continuum subtraction process but may still impact the 
inferred mole fractions indirectly though degeneracies in 
the surface gravity and effective temperature estimation. 
Uncertainties in T^ff translate into a small 0.1 dec) 
changes in the mole fractions. Future refinement of the 
mole fractions will benefit most by improved determina¬ 
tions of \og{g). Surface gravity broadens spectral lines 
and, in principle, the width of auto-correlation functions 
for models of various surface gravities could be used to 
determine gravity’s contribution to the total line broad¬ 
ening. Unfortunately, even a resolution of 4000 (^75 km 
s“^) is too low for such an exercise to yield better gravity 
estimates than those based on model fits to lower resolu¬ 
tion data (e.g., fitting shapes of H and K band spectra). 

4.5. Estimating K^z 

Vertical mixing is often characterized by the coeffi¬ 
cient for eddy diffusion (ATzz) where the vertical mixing 
timescale above the convection zone is Leff/and Lgff 
is an effective len gth scale usu ally a few tenths the pres¬ 
sure scale height (|SmithlIl9^ . Below the point in the 
atmosphere where the mixing timescales are shorter than 
chemical reaction timescales, the atmosphere will be in 


Figure 8. Top-. Photospheric mole fractions for methane as func¬ 
tions of the coefficient for eddy diffusion in models with \og{g) = 3.5 
and Tefj = 900, fOOO, and HOOK (solid, dashed and dotted lines, 
respectively). The horizontal bar indicates the range of methane 
mixing ratios inferred from the i^-band spectrum (Table 1). 

chemical equilibrium. Most chemical reactions are fast, 
however, the reactions governing CO and CH 4 have one 
or more ra te-limiting steps that result in long chemica l 
timescales (|Visscher et al.lIMT^ IZahnle fc MarlevI 1201411 
that rapidly increase with decreasing density, quickly ex¬ 
ceeding the age of the planet by many orders of mag¬ 
nitude. The ultimate consequence is that, for plausible 
values of K^z, the photospheric mole fractions of CO and 
CH 4 (as well as others, e.g., N 2 , NH 3 , and CO 2 ) may no 
longer depend on the photosphere conditions but instead 
on the conditions deeper in the atmosphere where the 
chemical and mixing timescales become comparable. 

The detection of CO in Jupiter’s atmosphere is an ex¬ 
cellent example of vertical mixing in action. For the very 
low temperatures in Jupiter’s atmosphere (~ lOOK), the 
majority of carbon is in CH 4 at an abundance that is 
essentially constant with height. Therefore, if quenching 
occurs in the atmosphere, the photospheric CH 4 will be 
unchanged. The equilibrium photospheric mole fractions 
of CO, on the other hand, are small and rapidly decreas¬ 
ing with height. If pure chemical equilibrium persisted 
throughout, CO would be nearly impossible to observe. 
Yet, photospheric CO absorption has been measured, ev¬ 
idence that non-equilibrium chemistry is at work, and 
the CO abundance used to estimate the value of K^z 
(|Feglev fc Lodder^ll994D . 

In the atmospheres of young giant exoplanets as hot as 
HR8799 b (^ lOOOA") the situation for CO and CH 4 i s 
reversed (j Barman et al.ll2nilbl : IZahnle fc MarlevI l2014fl . 
At pressures (between 1 and 10 bar, the temperatures are 
hot enough for CO to be thermochemically favored over 
CH 4 . Also, the quenching can occur at or near depths 
where equilibrium mole fractions are height-dependent 
for CH 4 and nearly height-independent for CO. An im¬ 
portant outcome is that, much like CO in Jupiter, mea¬ 
surements of CII 4 can be used estimate K^z in the atmo¬ 
spheres of planets like HR8799 b. 

The photospheric CII 4 mole fractions versus K^z are 
plotted in Figure | 8 ] for solar composition and a range of 
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Figure 9. Continuum-subtracted spectrum for HR8799 b (bla ck) compared to models (red) with low (top), high (bottom) and best¬ 
matching (middle) carbon-to-oxygen (C/0) ratios found in ^4.61 At low C/0, H 2 O and CO absorption lines are stronger than observed, 
while CH 4 absorption lines are weaker. At high C/0, CH 4 absorption is stronger than observed while H 2 O and CO are too weak. The 
short-wavelength portion of A-band is not shown for clarity, though the general level of agreement is comparable to what is plotted. 


Teff appropriate for HR8799 b. Only for \og{Kzz) > 10 
does K 22 become independent of CH 4 . The range of 
inferred mole fractions for methane (horizontal shaded 
region in Fig. 0 ) indicate \og{Kzz) around 7 for Teff= 
lOOOiF; however, the mole fractions at the quenching 
depth are temperature and gravity dependent leading to 
an uncertainty in of about 100 cm^ s“^. Large de¬ 
partures from solar C and O abundances could change 
the inferred K^z, however, as argued below, it is unlikely 
that C/0 is very far from solar. 

4.6. C/0 Ratio 

The C/0 ratio for HR8799 c was found to be slightly 
above that of the host star, tentatively favoring the core¬ 
accretion formation scenario over gravitational instabil¬ 
ity (K13). Determining the C/0 for HR8799 b is an 
important next step in understanding the formation his¬ 
tory of this planetary system. 

The photosphere of HR8799 b is not in chemical equi¬ 
librium, as discussed above, with both CO and H 2 O 
quenched deep in the atmospheres. Deep quenching of 
carbon and oxygen-bearing molecules has important im¬ 
plications for inferring the relative C and O element 
abundances from the molecular mole fractions. Even 
young giant planets, still hot from recent formation, have 
atmospheres cool enough to allow solids and liquids to 
form, including silicate grains (e.g., MgSiOa, Mg 2 Si 04 ). 
The abundance of these grains is limited by the overall 
abundance of Si, and will sequester a non-negligible frac¬ 


tion of the Oxygen atoms. The photosphere of HR8799 b 
is well below the temperatures needed for condensation 
to occur and, as a result, the inventory of Oxygen should 
account for both the mole fractions of silicate grains 
and Oxygen-bearing molecules. However, for the atmo¬ 
sphere of HR8799 b, models indicate that the tempera¬ 
ture at the quenching depths is above 2000K and, thus, 
above the condensation temperature for major oxygen- 
depleting grains. In the case of the HR8799 planets, 
non-equilibrium chemistry simplifies the C/0 ratio de¬ 
pendence on atmospheric mole fractions (N) to 

C N{CH4) + N{C0) 

O ~ N{H20) + N{CO) ’ ^ ’ 

and for small amounts of CH 4 , the C/0 ratio is deter¬ 
mined by H 2 O and CO alone. 

Equation 1 and the mole fractions determined above 
(Table 1) for each model atmosphere Teff and log(g), 
results in C/0 values between 0.4 and 0.7. This large 
range in C/0 is mostly driven by the comparably large 
range of acceptable values for log(g), with C/0 decreas¬ 
ing as log((/) increases. The mole fractions for the best 
matching model yield C/0 = 0.61 ±0.05 (with l-cr error 
determined from the formal errors in Table 1). Figure 
[5] compares the continuum-filtered observations to the 
best overall model from il4.6l (C/0 = 0.61) and models 
with high and low C/0, to illustrate the changes in the 
spectrum as the relative molecular abundances change. 
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Figure 10. Carbon and Oxygen abundances (solid lines) relative 
to stellar values, along with the corresponding C/O ratio (dashed 
line), using the lOberg et a DiEqU) model for planets forming be¬ 
tween the H 2 O and CO 2 frostlines. The disk is assumed to have 
the same grain/gas fraction as the interstellar medium (O.Of). See 
lOberg et al.l II2011I1 for the specific details of their model. 

5. ATMOSPHERIC COMPOSITION AND FORMATION 

If giant planets form primarily by a quick one-step 
process via gravitational instabilities (GI), their atmo¬ 
spheres s hould have element abun dances equal to the 
host star (iHelled fc Schubertll2009ll . On the other hand, 
if giant planets form primarily by the multi-step core¬ 
accretion (CA) process, a range of element abundances 
are possible (jOberg et al.ll20lill . The abundances of a 
gas giant’s atmosphere formed via CA primarily depend 
on the location of formation relative to the frostlines for 
major carbon and oxygen bearing molecules in the disk 
(namely H 2 O, CO 2 and CO) and the amount of solids ac¬ 
quired by the planet during the runaway accretion phase. 
The four planets orbiting HR8799 offer an excellent op¬ 
portunity to test this idea. Each planet currently orbits 
between the H 2 O and CO 2 frostlines and potentially built 
up atmospheres from gas with similar amounts of solids. 

In K13, the observed continuum-filtered spectrum of 
HR8799 c was compared to atmosphere models restricted 
to a sequence of C and O element abundances derived 
from the lOberg et al.l (|20Ilh chemical model. The Oberg 
et al. model provides values for the C and O abundances, 
relative to the stellar values, for different amounts of solid 
accretion during the buildup of the planetary envelope. 
These abundances are plotted in Figure [in] for planets 
forming between the H 2 O and CO 2 frostlines. In this 
model, both C and O abundances are linear functions of 
solid accretion, Msoiid/Mgas, with slope and intercept set 
by the fraction of C (or O) sequestered by condensate 
formation and the overall grain/gas fraction in the disk. 
These assumptions are based on observations of proto¬ 
planetary disks and the interstellar medium (see Table 
1 of Oberg et al. 2011). Solar C and O abund ances 
have been suggested for HR8799 (iSadakand 120061) an d 
are adopted as the baseline here (jAsolund et al.l 1200^ . 
It should be noted, however, that this star is a A-boo 
type star with solar C, N and O abundances but sub¬ 
solar Fe-peak elements. 

A high level of solid accretion during planet formation 



c/o 


Figure 11. Top-. Distribution of Ax'^ for models with C and 
O abundances following lOberg et all 11201111 compared to the full 
if-band spectrum with flux calibrated continuum for HR8799 b 
(solid line) and HR8799 c (dashed line). Bottom-. Same as top 
panel but for comparisons to the filtered (continuum-subtracted) 
if-band spectrum. The dotted line is the Ax'^ distribution for 
HR8799 c from K13 and the dashed line is the distribution for a 
reanalysis discussed in the text. 

raises both C and O abundances, with O abundances in¬ 
creasing more rapidly than C for the simple reason that 
H 2 0 -ice is the most abundant solid between the H 2 O and 
CO 2 frostlines. The combined effect is C/O decreasing 
as Mgoiid/Mgas increases. The model proposed by Oberg 
et al. is a simple prescription for a complex process and, 
consequently, deviations from this model are to be ex¬ 
pected. Despite its simplicity, the predicted C and O 
abundances provide an ideal baseline for testing poten¬ 
tial outcomes of CA formation specific to the HR8799 
system. 

In order to make a direct comparison to the K13 re¬ 
sults, the observed spectrum of HR8799 b was analyzed 
in a similar manner as HR8799 c. Given the CH 4 up¬ 
date made to the model atmospheres, the fit was re¬ 
peated for HR8799 c. Only the C/O values for atmo¬ 
sphere accretion occurring between the H 2 O and CO 2 
frostlines (the current locations of all four HR8799 plan¬ 
ets) were used (see Fig. (TO]). A was calculated for each 
continuum-subtracted synthetic spectrum in the and 
log(( 7 ) range described above. The resulting distribu¬ 
tions are plotted in Figure ITT] The best-matching C/O 
for HR8799 b is 0.66/)Q gg and the revised best-matching 
C/O for HR8799 c is 0.64 /)oud closely matching K13, 
but with a broader distribution. The new ExoMol 
methane absorption stre ngths across the AT-b and are 
lower than those from the lWarmbier et'aT] ([200911 linelist 
used in the KI3 analysis resulting in smaller values 
for larger C/O and, hence, a broader distribution of ~x^. 

The AT-band continuum shape can also be used to esti¬ 
mate the relative abundance of H 2 O and CO. Given the 
wider angular separation of b than c (1.7 versus 1") from 
the star, the continuum of b is less affected by residual 
speckles. The observed spectrum, with continuum intact, 
was compared to the same set of atmosphere models and 





































10 


Barman et al. 



Figure 12. Photospheric mole fractions of CO, H 2 O, and CH 4 
as functions of C/0 for atm osphere models with C and O abun¬ 
dances as plotted in Fig. 1101 For each molecule and C/O, a range 
of model values is possible given the range of effective temperature 
and gravity explored here (900 - HOOK, 3.0 - 4.0). The ranges 
of mole fractions are plotted as shaded regions for H 2 O and CH 4 . 
The range of CO mole fractions is narrow and plotted as a black 
line with varying thickness. These regions correspond to the photo- 
spheric molecular abundances specifically for models of HR8799 b, 
and the assumptions made here about vertical mixing of CO and 
CH 4 . The best-matching model with Oberg et al. based abun¬ 
dances has a C/O ~ 0.66 (see ^ and indicated by the vertical 
dotted line. The full range of allowed C/O from not shown 

to avoid cluttering the Figure, but extends from 0.54 to 0.7. The 
intersections of this dotted line and the shaded regions indicate 
the molecular mole fractions at the photosphere of the best-fitting 
model. For comparison, solar abundances and C/O are indicated 
by the vertical dashed line. T he m ole fractions and best matching 
range of C/O determined in ^4.61 are plotted as red symbols with 
horizontal and vertical error-bars. The different model fitting pro¬ 
cedures result in consistent mole fractions for H 2 O and CH 4 , while 
the scaling procedure results in CO lower than in all of the model 
atmospheres with Oberg et al. based C and O abundances. 

C and O abundances. The distribution from this com¬ 
parison is plotted in the top panel of Fig. [TTland shows 
both planets having C/O between ~ 0.55 and 0.7, for 
similar Ax^ as in the continuum filtered comparison. 

The mole fractions of H 2 O, CO and CH 4 in the pho¬ 
tosphere of the best-matching CA-specific atmosphere 
model can be compared to those found in H4.41 where 
the molecular abundances were determined independent 
of any formation model. The mole fractions for the CA- 
specific atmosphere models are plotted in Figure [T^ as 
functions of C/O, for the range of Tes and log{g) ex¬ 
plored here. Solar (non-equilibrium) mole fractions and 
solar C/O are indicated in this Figure by the vertical 
dashed line. The CA-specific atmosphere model with 
C/O ^ 0.66 has H 2 O and CH 4 mole fractions that closely 
match those found in >14.41 These values are plotted as 
red symbols with error bars estimated from the formal 
uncertainties in Table 1 and the CO abundances cor¬ 
responding to the 1—( T ran ge in C/O. CO, however, is 
higher than found in >14.41 for all C/O explored in the 
CA-specific atmospheres. Even if log( 5 ) and T^ff are var¬ 
ied (see Table 1), the CO found by freely scaling the 
abundances remains below what is predicted in the CA- 
specific atmospheres. Lower K^z would lower CO, but 


doing so would quickly increase CH 4 to values well above 
what is observed. A lower CO could indicate additional 
chemical processes not accounted for in the atmosphere 
or formation models. The significance of this CO differ¬ 
ence is hard to assess given the various source of uncer¬ 
tainty in the models not accounted for in the formal mole 
fraction error-bars. It is possible, for example, that the 
CO line broadening in the atmosphere models is under 
(or over) estimated, resulting in systematic errors in the 
inferred abundances (this will be investigated in a future 
paper). 

Compared to H 2 O and CH 4 , CO has the weakest de¬ 
pendence on C/O. CH 4 and H 2 O, however, are fairly 
sensitive to C/O, but their mole fractions can be small 
(< 10“'*) for super-solar C/O. Measuring C/O from CO 
alone would require a CO abundance more precisely de¬ 
rived than most data presently allow. Perhaps of greater 
interest, however, are the ratios of the various molecules 
as C/O changes. For example, for all but the smallest 
C/O, CO/H 2 O should be greater than 1. Only for low 
surface gravity (< 3.5) do the model fits from > 14.41 re¬ 
sult in N(CO) > N(H 20 ) and, as gravity increase, the 
discrepancy between the predicted and observed CO and 
H 2 O abundances grows. This may provide indirect evi¬ 
dence supporting low surface gravity for HR8799 b. 

6 . SUMMARY AND CONCLUSIONS 

A new AT-band spectrum of HR8799 b, with spectral 
resolution of ^ 4000, was measured and simultaneous 
detections of water, carbon monoxide and methane ab¬ 
sorption lines were made. Identification of lines from 
each molecule is possible by eye and confirmed by 
cross-correlating the observed spectrum with absorption 
templates for each molecule. The abundance of each 
molecule was determined by a fitting procedure where 
the mole fractions are treated as free parameters in se¬ 
quences of synthetic spectra. The ensemble of near-IR 
to thermal IR observations for HR8799 b are best re¬ 
produced by an atmosphere with Teg= 1000 , log( 5 ) = 
3.5. The log mole fractions are found to be between - 
3.09 and -3.91 for H 2 O, between -3.30 and -3.72 for CO, 
and between -5.06 and -5.85 for CH 4 . The best match¬ 
ing models have C/O between 0.55 and 0.7 depending 
heavily on log((;), as also found in K13, with C/O de¬ 
creasing with increasing log(( 7 ). If the surface gravity 
is closer to log(< 7 ) = 4, as pre dicted by hot-start cool¬ 
ing tracks (jBaraffe et al.ll200'll . then CO would be less 
abundant than H 2 O (see Table 1) with a C/O ~ 0.4. 
Such a situation would require C and O abundances that 
deviate from those predicted by the lOberg et al.l (|201iri 
model, perhaps resulting from planetesimal accretion or 
core dredging. 

A cross-correlation analysis of i7-band data from Bll, 
taken at a similar spectral resolution as the X-band data, 
revealed no methane signature. A detection of CH 4 ab¬ 
sorption lines in the X-band but not in H may seem un¬ 
expected when compared to late-type brown dwarfs that 
often show more prominent CH 4 absorption in H than 
X. As shown here, the methane abundance in HR8799 b 
is many orders of magnitude below that found in the 
atmospheres of field brown dwarfs, resulting in weaker 
methane absorption at all wavelengths. For photospheric 
temperatures and pressures appropriate for HR8799 b, 
the average methane opacity (in units of cm^/molecule) 
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across K is about 10 times stronger than the average 
opacity across H, while the average water opacity is 
about 10 times stronger than methane in H and 10 times 
weaker than methane in K. Consequently, water opacity 
likely overwhelms that of methane across the H band. 
The characteristics of these opacity sources, combined 
with lower SNR in the i7-band data, all contribute to 
the non-detection of methane in the i7-band spectrum. 

Low methane abundance in HR8799 b is consistent 
with previous studies of non-equilibrium chemistry in 
hot, low gravity atmospheres. The mixing ratio of CH 4 
in young giant exoplanets is, potentially, a useful probe of 
the vertical mixing. Using the inferred CH 4 abundance, 
the coefficient of eddy diffusion was found to be greater 
than 10 ® and very near 10 ^ cm^/s for the preferred Teff = 
lOOOK. These ATzz values are consistent with values ap¬ 
propr i ate for Jupiter’s d e ep tro posphere (IVisscher et al.l 
1201011 . IZahnle fc MarlevI (|2014ll explore the topic of non¬ 
equilibrium CH 4 in the context of young giant planets 
and their predictions are in good agreement with the 
values presented here. 

HR8799 b and c are the first exoplanets to have C/0 
ratios determined from spectroscopic data with high SNR 
and high spectral resolution (i? ^ 4000). Young giant 
planets, still orbiting their star beyond the H 20 -frostline, 
provide important opportunities to study the link be¬ 
tween formation and atmospheric chemistry. The C/0 
ratios of both HR8799 b and c are similar and poten¬ 
tially super-stellar, but stellar values are not completely 
excluded. Refinements of C/0 and individual molecu¬ 
lar abundances will require improvements in the deter¬ 
mination of surface gravity, independent of evolution¬ 
ary models, perhaps with near-IR observations at even 
higher spectral resolution. Observations with the Hub¬ 
ble Space Telescope, at wavelengths within near-IR water 
absorption bands, might also help narrow the range of 
allowed abundances found here (Rajan et ah, in prep). 
But comparisons to such data would face similar limi¬ 
tations as here, associate with effective temperature and 
gravity. More observations across the thermal infrared (3 
to 5/rm), perhaps with future ground-based integral field 
spectrographs or with the James Webb Space Telescope, 
could provide more complete coverage of strong CH 4 , CO 
and CO 2 absorption bands. New measurements across 
this wavelength range would also help improve effective 
temperature and gravity estimates. Model uncertainties 
need to be better understood as well, in particular those 
associated with clouds and natural molecular line broad¬ 
ening. 

There are now three young directly imaged planet- 
mass companions (HR8799 b, c and 2M1207B) that show 
evidence of disequilibrium chemistry in their deep tropo¬ 
sphere. While three is too few to constitute a trend, 
the data and models suggest that chemical quenching is 
common among this class of exoplanet. For planets as 
warm as those orbiting HR8799, the quenching likely oc¬ 
curs below the condensation depths and, consequently, 
bulk C/0 ratios are inferable form their photospheric 
H 2 O and CO abundances. Given the broad absorption 
features of these molecules, it may be possible to infer 
C/0 for massive planets from the low resolution spectra 
provided by CPI, P1640, and SPHERE. 
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